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The origin of the antiferromagnetic order and puzzling properties of LaTiOa as well as the 
magnetic phase diagram of the perovskite titanates are studied theoretically. We show that in 
LaTiOs, the t2g degeneracy is eventually lifted by the La cations in the GdFeOa-type structure, 
which generates a crystal field with nearly trigonal symmetry. This allows the description of the 
low-energy structure of LaTiOs by a single-band Hubbard model as a good starting point. The 
lowest-orbital occupation in this crystal field stabilizes the AFM(G) state, and well explains the 
spin- wave spectrum of LaTiOa obtained by the neutron scattering experiment. The orbital-spin 
structures for .RTiOs with R—Pr, Nd and Sm are also accounted for by the same mechanism. 
We point out that through generating the R crystal field, the GdFeOa-type distortion has 
a universal relevance in determining the orbital-spin structure of the perovskite compounds 
in competition with the Jahn- Teller mechanism, which has been overlooked in the literature. 
Since the GdFeOa-type distortion is a universal phenomenon as is seen in a large number of 
perovskite compounds, this mechanism may also play important roles in other compounds of 
this type. 
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1. Introduction 

A series of the perovskite-type transition-metal oxides 
have attracted much interest since strongly correlated 
electrons in some of these compounds show astonishing 
transport phenomena upon carrier doping such as high Tc 
superconductivity in layered cuprates and colossal mag- 
netoresistance in manganites.^^ Perovskite titanium ox- 
ide is also one of the typical examples. However, in spite 
of the structural similarity, the titanates show very dif- 
ferent properties from the above two compounds. 

The perovskite titanate Ri-xAxTiOs (i?=rare-earth 
ions, A=alkaline ions) shows a transition from insulator 
to metallic state upon carrier doping achieved by increas- 
ing a;.^"^^ This metallic phase proved to be well described 
by the Fermi liquid picture, and the origin of the very 
different behavior of the titanates from the cuprates and 
manganites has been a subject of very intensive studies. 
More concretely, strong mass enhancement was found in 
the specific heat and the susceptibility near the metal- 
insulator transition point, and discussed in connec- 
tion to the nature of the Mott transition. ^"^^^ The nature 
of the transition shows a marked contrast with the cases 
in the cuprates and manganites, and the titanates have 
been recognized as a touchstone material in clarifying the 
strong electron correlation effects. 

In order to clarify the nature of the metallic phase and 
the transition, it is also essentially important to clarify 
the electronic structure of the end Mott-Hubbard insu- 
lator i?Ti03. In fact, i^TiOa shows a number of puz- 
zling and rich electronic structures arising from the or- 
bital degrees of freedom. In these compounds, Ti''+ has 
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Fig. 1. Crystal structure of -RTiOa with the GdFeOs-type distor- 
tion. 

a configuration in which one of the threefold t2g or- 
bitals is occupied by an electron, and the 3d t2g bands 
are at the Fermi level in contrast with the cuprates and 
manganites with eg bands at the Fermi level. One of the 
central puzzles in i?Ti03 is, as we will describe details be- 
low, that orbitally degenerate systems, in general, have a 
strong tendency towards antiferro-orbital and ferromag- 
netic state, ^' while LaTiOs shows actually antiferromag- 
netic order at low temperatures. Here, "antiferro-orbital 
order" is defined as an order of orbital polarization with 
the staggered pattern in analogy with the staggered spin 
polarization in the antiferromagnetic order. One of the 
purposes of this paper is to solve this puzzle and clarify 
the origin of the antiferromagnetic order. In the course of 
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the clarification, we show that LaTiOa and related com- 
pounds are, in reality, a good prototype of the single- 
band Hubbard model on a cubic lattice. 

The crystal structure of i?Ti03 is a pseudocubic per- 
ovskite with an orthorhombic distortion (GdFeOa-type 
distortion) in which the TiOe octahedra forming the per- 
ovskite lattice tilt alternatingly. Note that, as a first ap- 
proximation, the cubic TiOg octahedra are not atom- 
ically distorted with this distortion. In this structure, 
the unit cell contains four TiOg octahedra (sites 1-4) as 
shown in Fig. 1. The magnitude of the distortion depends 
on the ionic radii of the R ions. With a small ionic radius 
of the R ion, the lattice structure is more distorted and 
the Ti-O-Ti bond angle is decreased more significantly 
from 180°. For example, in LaTiOa, the bond angle is 
157° (a5-plane) and 156° (c-axis), but 144° (a6-plane) 
and 140° (c-axis) in YTiOs.") The ionic radii of La and 
Y ions are 1.17 A and 1.04 A, respectively. Here, we note 
that the magnitude of the distortion can also be con- 
trolled by using solid-solution systems. For example, by 
varying the Y concentration in Lai_a;Yj.Ti03, we can 
control the bond angle almost continuously from 157° 
{x=0) to 140° (a;=l). In the perovskite titanates, this 
bond angle distortion controls the interplay of the or- 
bital, spin and lattice degrees of freedom, and the system 
shows various orbital-spin phases and their phase transi- 
tions due to their interplay. It is widely recognized that 
the bond angle primarily controls the electron transfers 
between the neighboring Ti t2g orbitals mediated by the 
O 2p orbitals, thereby reduces the t2g bandwidth with 
decreasing bond angle. In this paper, we will show that 
the tilting plays another crucial role in determining the 
electronic structures. 

Recently, a magnetic phase diagram in the plane of 
temperature and the magnitude of this distortion was 
obtained experimentally for iiTiOa, which exhibits an 
antifcrromagnctic (AFM)-to-fcrromagnctic (FM) phase 
transition (see Fig. 2 upper panel). LaTiOa with 
the smallest distortion shows a G-type AFM (AFM(G)) 
ground state, in which spins arc aligned antifcrromagnct- 
ically in all y and z directions as shown in the inset 
of the upper panel of Fig. 2. The magnetic moment is 
0.45 /J-B:^^' which is reduced from spin-1/2 moment. The 
Neel temperature (Tn) is about 130 K. With increas- 
ing GdFeOs-type distortion, Tn decreases and is strongly 
depressed at SmTiOa, subsequently a FM ordering ap- 
pears. In the significantly distorted compounds such as 
GdTiOs and YTiOs, a FM ground state accompanied by 
a large .Jahn- Teller distortion is realized. A similar phase 
diagram was also obtained for Lai_a;Yj,Ti03 (see Fig. 2 
lower panel). ^^'^^^ 

A suppression of the Curie temperature (Tq) around 
the AFM-FM phase boundary is understood from the 
strong two-dimensional anisotropy in the FM coupling 
near the transition point. '^^^ On the other hand, in 
spite of a lot of attempt and effort, the origin of the 
AFM(G) structure in LaTiOs has been controversial. We 
briefly summarize the history of the study on this con- 
troversy in the following. 

Previous diffraction studies^^' showed that a Jahn- 
Teller type distortion of the TiOe octahedra in LaTiOs 
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Fig. 2. Experimentally obtained magnetic phase diagram for 
-RTiOa (upper panel, Ref. 12, 13), and that for Lai_a:Y2,Ti03 
(lower panel, Ref. 15, 16). In these phase diagrams, Tn and Tq 
are plotted as functions of the ionic radius of the R ion and 
the unit cell volume, respectively. (The lines arc guides for the 
eyes.) Note that the unit cell volume of Lai—j, YajTiOa is ap- 
proximately proportional to the Y concentration (x), and well 
characterizes the magnitude of the GdFcOs-type distortion. 
Here, the closed symbols and the open symbols indicate Tn and 
Tq , respectively. The values of the unit cell volume are deduced 
from the data obtained by the previous x-ray diffraction mea- 
surement (Ref. 15). 



was undetectablly small in contrast with the cuprates 
and manganites. If this is true, in LaTiOa, the crystal 
field from the O ions surrounding a Ti^"*" ion has a cubic 
symmetry so that the degeneracy of the t2g orbitals is 
expected to survive. In fact, a number of papers have 
assumed this degeneracy as we describe below, and this 
is a source of controversies. 

Under this circumstance, the AFM(G) ground state 
has been surprising and controversial since in the or- 
bitally degenerate system, it is theoretically expected 
that a FM state with antiferro-orbital ordering is fa- 
vored both by electron transfers and by Hund's rule cou- 
pling. Indeed, a recent theoretical study based on a 
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Kugcl-Khomskii model showed that in the cubic lattice 
with threefold degenerate t2g orbitals, a FM state is al- 
most always stable, and the AFM(G) state appears only 
in the unphysical regions of spin-exchange parameters. 

Mizokawa and Fujimori performed a Hartree-Fock 
analysis of the multiband d-p model, and proposed that 
the spin-orbit (LS) interaction lifts the t2g degeneracy, 
and the AFM(G) state in LaTiOs is accompanied by 
the LS ground state. ^^'^^^ In the LS ground state, un- 
quenched orbital moment antiparallcl to the spin mo- 
ment is necessarily induced. At first sight, this seems to 
be consistent with the experimentally observed reduced 
magnetic moment in this compound. 

However, in contrast with this naive expectation, a 
recent neutron scattering study has shown a spin-wave 
spectrum well described by an isotropic spin- 1/2 Heisen- 
berg model (J ~15.5 meV) with a considerably small 
spin gap.^^^ This indicates that the orbital moment is 
almost quenched, and the LS interaction is irrelevant 
since if the orbital moment exists, an induced Ising-type 
anisotropy in the spin sector has to generate a large 
spin gap. Actually, an exact diagonalization study has re- 
cently shown that the Kugel-Khomskii model with LS in- 
teraction does not describe the observed spin-wave spec- 
trum. In addition, a model Hartrce-Fock study done 
by one of the authors^^^ also showed that the LS ground 
state cannot reproduce the AFM(G) ground state by ex- 
amination of a solution which was overlooked in the pre- 
vious Mizokawa and Fujimori's study. In Ref. 31, it was 
shown that even without the static Jahn- Teller distor- 
tion, a FM state out of which two states {yz + izx) t 
and xy t are alternating is stabilized both by the LS 
interaction and by the spin-orbital superexchange inter- 
action. 

To explain the neutron scattering result, a possible 
orbital liquid state was proposed on the basis of small 
orbital-exchange interaction in the AFM(G) spin struc- 
ture. ^^"^'^^ In these theories, the AFM(G) ordering is as- 
sumed a priori. However, in the titanates, the spins and 
orbitals strongly couple to each other, and both degrees 
of freedom cannot be determined independently. There- 
fore, the origin of the AFM(G) state in LaTiOs is to be 
clarified in a self-consistent manner. More importantly, 
a FM state with antiferro-orbital ordering was theoreti- 
cally expected to be more stable in this system. Thus, the 
assumption of the AFM(G) order on this basis is hard 
to be justified since the relative stability of the AFM(G) 
to FM state has never been examined. In addition, a 
recent heat capacity measurement showed that the most 
of the low-temperature heat capacity arises from magnon 
contributions, which contradicts the prediction from the 
orbital liquid theory about the orbital contributions. 

A possible trigonal {D34) distortion of the TiOg octa- 
hedra was proposed to explain the emergence of AFM(G) 
structure. With this distortion, the TiOe octahedron 
is contracted along the threefold direction, and the three- 
fold degenerate t2g levels split into a nondegenerate lower 
aig level and twofold-degenerate higher Cg levels (see 
Fig. 3). Occupations of the aig orbitals well explain the 
emergence of the AFM(G) ordering and the isotropic 
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trigonal 

Fig. 3. With the D,j^ distortion, the TiOg octahedron is con- 
tracted along a trigonal direction. There exist two kinds of O-O 
bonds, shorter O-O bond and longer O-O bond as presented 
by solid lines and dashed lines, respectively. As a result of this 
distortion, the threefold degenerate t2g levels split into a nonde- 
generate lower aig level and twofold degenerate eg levels. With 
the D-^i^ crystal field with [1,1,1] trigonal axis, the representation 
of the aig orbital is -y=(xy + yz + zx). 



spin-wave spectrum. However, this distortion had not 
been observed clearly upon theoretical proposal. 

NMR studies have also applied to this compound. Pre- 
viously, the NMR spectrum was analyzed based on the 
orbital liquid picture. However, it has recently been 
claimed that there exists a discrepancy. Indeed, it was 
recently clarified that the spectrum is well described by 
an orbital ordering model proposed in the above D^d- 
distortion scenario instead of the orbital liquid model. 
Moreover, a recent resonant x-ray scattering result also 
indicates the orbital ordering in the series of the AFM(G) 
compounds of i?Ti03 {R =La, Pr, Nd and Sm), which 
has the same symmetry with that observed in the ferro- 
magnetic YTiOa and GdTiOs,^®^ which also contradicts 
the orbital liquid picture. 

Several first-principles methods were also applied to 
LaTiOs. However, LDA and GGA failed to reproduce 
not only AFM(G) structure but also insulating behav- 
ior since the strong electron correlation is insufficiently 
treated. In addition, although the insulating gap 
was reproduced, the AFM(G) structure was not obtained 
in LDA-|-i7 calculation."'^' These indicate that in this 
compound, the effect of the strong electron correlations 
is crucially important, which is of general interest of the 
condensed matter physics. 

In such ways as discussed above, the problem of the 
orbital-spin structure in LaTiOs has been studied inten- 
sively, but a consistent theory has not emerged, and is 
still under hot debates. In fact, when we were perform- 
ing the calculations and were preparing this manuscript, 
several reports of new experimental findings and theo- 
retical results appeared in the preprint server, which we 
will introduce and discuss later in the context. 
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To fully understand electron correlation effects in the 
perovskite compounds and in more general, in d-electron 
systems, it is recognized that puzzling properties of 
the Ti perovskites have to be explained by a convinc- 
ing theory. The Ti perovskite compound also offers 
touchstone materials for the understanding of orbital 
physics, which is a hot topic of the physics of electron 
correlation."^^) Because of the threefold degeneracy of d} 
state in the titanates, it is expected that the complete 
understanding of the Ti perovskites contributes in 
clarifying roles of orbital degrees of freedom as well 
as coupling and interplay between magnetism and 
orbitals in the strongly correlated electron systems. The 
purpose of the present theoretical study is to solve the 
long-standing puzzles in the perovskite titanates and to 
offer an overall understanding of the interplay in these 
compounds. 

In the perovskite compounds, the Jahn- Teller distor- 
tion often plays crucial roles in determining low-energy 
electronic state by lifting the orbital degeneracy. How- 
ever, with early transition-metal ions, the .Jahn- Teller 
coupling is considerably weaker than that in the late- 
3d compounds. Indeed, in contrast with the cuprates and 
manganites, a sizable Jahn- Teller distortion had not been 
observed in LaTiOa as already mentioned above. With- 
out the Jahn- Teller distortion, the crystal field from the 
ligand O ions acting on the Ti^+ ion has a cubic symme- 
try. 

In addition to the Jahn- Teller distortion, the GdFeOs- 
type lattice distortion is another generic phenomenon in 
perovskites. However, it has been assumed mainly to con- 
trol the bandwidth through the M-O-M angle variation 
with M being a transition-metal ion,^^''*'^) while its di- 
rect effects on interplay of spins and orbitals have not 
been considered seriously. In fact, the knowledge about 
the symmetry lowering caiiscd by the GdFeO.j-typc dis- 
tortion in principle could lead to the idea for the lift of 
the f2g-orbital degeneracy. However, in contrast with the 
Jahn- Teller distortion, a imit TiOg octahedron as well 
as its cubic symmetry are not distorted in the generic 
GdFeOa-typc distortion so that it has been assumed in 
most of the literature that the threefold degeneracy of 
the t2g orbitals is retained even with this distortion. Ac- 
tually, the recent orbital physics is usually based on the 
degenerate t2g and Cg orbitals unless the Jahn- Teller type 
octahedral distortion exists. However, this assumption is 
in fact a source of the puzzles recognized in LaTiOs. 

To solve the puzzles and controversies in the titanates, 
we consider that it is important to reveal how the basic 
degeneracy of the orbitals is retained or split in the real- 
istic experimental situation since this is a starting point 
of studies on orbital physics for perovskite compounds. 
We believe that it is crucial to establish the precise es- 
timate about the degree of orbital degeneracy under the 
existence of the GdFeOa-type distortion while its effects 
have never been considered seriously. 

In this paper, we clarify that the generic GdFeOs-type 
distortion actually generates a new mechanism for con- 
trol of low-energy orbital-spin structure through lifting 
of the orbital degeneracy by generating crystal fields of 
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the R ions. This mechanism competes with that of the 
Jahn- Teller distortion and the LS interaction. 

When the Jahn- Teller distortion of the TiOe octahe- 
dra is small, we expect that the R crystal field plays 
a dominant role on the electronic structure in the per- 
ovskites by lifting the t2g degeneracy. By utilizing the 
positional parameters of the La ions measured by the x- 
ray diffraction study, we construct the Hamiltonian for 
the crystal field from the La ions in LaTiOs. Analyses of 
the obtained Hamiltonian show that the GdFeOa-type 
distortion, and resulting displacements of the La ions 
turn out to generate a crystal field with nearly trigonal 
symmetry. As a result, the threefold degenerate cubic-f2g 
levels split into three nondegenerate levels. Further, we 
study the stability and properties of the AFM(G) state 
in this crystal field by using an effective spin-pseudospin 
Hamiltonian constructed through the second-order per- 
turbational expansion in terms of the electron transfers 
in the limit of the strong Coulomb repulsion. The calcu- 
lations of energies and spin-exchange constant based on 
this effective Hamiltonian show that lifting of the t2g de- 
generacy and the lowest-orbital occupation well explain 
the emergence and properties of the AFM(G) state in 
LaTiOg. 

We also discuss a recent experimental finding of the 
octahedral distortion in LaTiOg. As already mentioned 
above, we previously predicted the D^d distortion of 
the TiOe octahedra in LaTiOg as an origin of the 
AFM(G) ordering. After this theoretical proposal, sev- 
eral diffraction studies have been done,*^"^'^' and re- 
cently, the distortion was actually detected by Cwik. et 
al. (Ref. 48). We discuss that the emergence of this dis- 
tortion is naturally explained by our theory, where the 
crystal field from the displaced R ions drives the distor- 
tion of the TiOe octahedron. We also show that in spite 
of this distortion, our approach in which the crystal field 
from the R ions is treated as a perturbation to the t2g 
degeneracy is justified. 

Further, in order to clarify the relation between the 
GclFeO.j-type distortion and the crystal field of the R 
ions, we also study the crystal field Hamiltonians for 
RTiOs with R being Pr, Nd and Sm. Analyses of the 
crystal field Hamiltonians again reproduce the AFM(G) 
state in each compound. The calculated spin-exchange 
constants well explain the observed reduction of Tn as 
the GdFeOg-type distortion increases. Results of a recent 
resonant x-ray scattering experiment are also discussed 
by analyzing the orbital structures in detail. In addition, 
wc also show that in SmTiOg, the R crystal field caused 
by the GdFeOg-type distortion strongly competes with 
the O crystal field due to the Jahn- Teller type TiOe dis- 
tortion. The orbital-spin structure and the depression of 
Tn in SmTiOg are discussed by considering this compe- 
tition. 

We point out the importance of the GdFeOg-type dis- 
tortion as a universal control mechanism of orbital-spin 
structure through generating the R crystal field in the 
perovskite compounds, which has been so far overlooked. 
This mechanism is substantially different from a usual 
Jahn- Teller mechanism. 

The organization of this paper is as follows. In Sec. 2, 
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we construct the effective spin-pseudospin Hamiltonian 
of the perovskite Ti oxides by following an approach sim- 
ilar to the Kugel-Khomskii formulation.^^' In Sec. 3, 
we examine the crystal field Hamiltonian of the La ions in 
LaTiOs. We discuss the origin and nature of the AFM(G) 
state as well as several experimental results for LaTiO.-j. 
In particular, wc focus on the spin-wave spectrum ob- 
tained by a neutron scattering study in detail. In Sec. 4, 
we examine the crystal field of the R cations in iJTiOs 
with R being Pr, Nd, and Sm. On the basis of this analy- 
sis, we discuss the behavior of Tn in the magnetic phase 
diagram and recent results of a resonant x-ray scatter- 
ing. In Sec. 5, wc discuss the following five points upon 
our proposal: (1) Available experimental results for the 
titanates are well explained within our theory. (2) The 
mechanism proposed in this paper is novel and univer- 
sal, which controls the orbital-spin structure in the per- 
ovskites, and is completely different from the usual Jahn- 
Tcllcr one. (3) Numerical results here obtained are suf- 
ficiently accurate even though they are calculated based 
on a point charge model with the nearest neighbor ions. 
Further, the predicted isotropy of the AFM(G) coupling 
as well as the properties of the orbital-spin states in 
RTiOs are robust. (4) The experimentally obtained mag- 
netic phase diagram is well explained in the light of the 
present results together with the previous ones. (5) The 
experimentally observed reduction of the magnetic mo- 
ment in LaTiOs is quantitatively understood when we 
consider a large amount of itinerant fluctuations in this 
compound. Section 6 is devoted to the summary. A short 
report on a part of the content of this paper has already 
been published. In this paper, we describe a compre- 
hensive framework and results on this issue together with 
additional new results and more detailed discussions. 

2. Microscopic Derivation of the Effective Spin- 
Pseudospin Hamiltonian 

In this section, we derive the effective spin-pseudospin 
Hamiltonian to describe the realistic system of the per- 
ovskite titanates in the limit of the strong Coulomb re- 
pulsion. We start with the multiband d-p model in which 
the full degeneracies of the Ti 3d and O 2p orbitals as 
well as the on-site Coulomb and exchange interactions 
are taken into account. The Hamiltonian is given by 

iJ'^P = HdO + Hp + Htdp + Htpp + i?on-site 



with 



(2.1) 



(2.2) 



(2.3) 



(2.4) 



Htpp= J2 t'jlj'i'P]iaPj'i''^ + ^-c-' (2-5) 

(2.6) 



where dl^^ is a creation operator of an electron with spin 

a{='\, l) in the 3d orbital 7 at Ti site i, and pj^^ is a 
creation operator of an electron with spin cr(=t,i) in 
the 2p orbital / at oxygen site j. Here, we choose the xy, 
yz, zx, 3z^ — and .x^ — orbitals as the basis of the 
3d orbitals. Hdo and Hp express the bare level energies 
of the Ti 3d and O 2p orbitals, respectively. Htdp and 
Htpp express the d-p and p-p transfers, respectively, ff^ ji 
and j/;/ are the nearest-neighbor d-p and p-p transfers 
given in terms of Slater-Koster parameters V^dx, Vpda, 
Vppjr and Vppa-^^'^^"^ The term iJon-sitc represents the 
on-site d-d Coulomb interactions, which consists of the 
following four contributions: 



(2.7) 



Hu' 



Hj 



E 



I'dl 



uai^^di^^dly^,diY„-, (2.8) 



E 



zja'^iYo-d^Y^idi^a-' J (2-9) 



Hf = E j'^Udiy^dl-yidii'l, (2-10) 

where and Hu' are the intra- and inter-orbital 

Coulomb interactions, respectively, and Hj and Hji de- 
note the exchange interactions. The term Hj is the ori- 
gin of the Hund's rule coupling. The term Hj/ gives the 
ti-pair hopping between the 3d orbitals on the same Ti 
atom. These interactions are expressed using Kanamori 
parameters, u, u', j and f which satisfy the following 
relations:^**' u = U + ^j, u' = u - 2j and j = j'. 
Here, U gives the magnitude of the multiplet-averaged 
d-d Coulomb interaction. The charge-transfer energy A, 
which describes the energy difference between the occu- 
pied O 2p level and the unoccupied Ti 3d level, is de- 
fined by using U and the energies of the bare Ti 3d and 
O 2p orbitals £° and Sp as A — e'^ + U — Sp since the 
characteristic unoccupied 3d level energy at the singly- 
occupied Ti site is Sd+U. The values of A, U and Vpda 
are estimated by the analyses of the LDA band calcu- 
lation^^) and the cluster-model analyses of the valence- 
band and transition-metal 2p core-level photoemission 
spectra. '''^"''^) In the present d-p Hamiltonian, wc have 
ignored the Coulomb interaction between two electrons 
on the O 2p orbitals as well as those on the O 2p and the 
Ti 3d orbitals. We take the values of these parameters as 
A = 5.5 eV, U = 4.0 eV, Vpda = -2.4 eV, Vpd^ = 1.33 
eV, Vppa- = 0.52 eV, Vpp^ = -0.15 eV and j = 0.46 eV. 

In the path-integral formalism, the expression of the 
partition function is given by 



^ = j 'Ddi^„'Ddi^„'Dpji„'Dpji^ exp 
with 



/3 



dTL(T) 



(2.11) 



= H,i -\- H,,i -\- Hi -\- H-ji , 



L{t) = H'^p{t) + J2 - ^)^' 



17(T 
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hi," 

where r denotes the imaginary time introduced in the 
path-integral formaUsm, and di-^a, di^a, pjia and pjia 
are the Grassman variables corresponding to the oper- 
ators dJ^CT' ^ii<^^ Pjia Pji<y^ respectively. By using the 
Matsubara- frequency representation; 

di^^ir) = ^ ^di^^K)e-''^"^, (2.13) 
pMr) = ^ ^p,^.(u;„)e-•'^"^ (2.14) 



we have 
Z 



dija{i^n)T^di^a{uJ„)'Dpji^^^^^)'Dpji^{u>n) 



X exp 



(2.15) 



with 



L{0Jn) = ^ di^a{-'^l^v 



ll)dr 



X] P3l<^i~^^n + £p - IJ')PjlcT 



3,1,'' 



i,'y.j,l.t7 

+ t^ljU'PJlTPj' I' a + C.C. 

3,l,3',l'cr 



(2.16) 



After integrating over the p and p, the partition function 
is rewritten as 



Z= 'Ddi^a{'^n)1^di^a{'^n)e:iq> 



where 



(2.17) 

(2.18) 



with 



- Z] dija-{ — iu)n + £d — ^ijdija 

+-f^on-site, (2-19) 
Ld2{c0n) = J2 ^i-y^ (2-20) 

Here, the matrix Hjij'ii{icOn) takes the form 
Hjij'i'iiuJn) = -{-ii^n + isp- IJ-Mi-j't' -t^ijn^- (2.22) 

Since the energy scale of electron transfers tpd (^1.3 
eV) and tpp (~0.3 eV) are relatively smaller than the 
characteristic energy of the on-site interaction U (=4.0 
eV), it may be justified to treat the term Ld2{'^n) as 
a perturbation to the iJon-site term. Therefore, we ap- 



proximate iujn in Ld2{^n) with poles of Ldi(w„), in which 
dynamics of the O 2p orbitals are integrated out. As a 
result, we obtain the expression of the effective super- 
transfer matrix clement between the Ti 3d orbitals, and 
the expression of the crystal field from the ligand oxygens 
acting on the Ti 3d orbitals as follows; 

= E H-iy'A4 + U-„)t'^.it%,.,,, (2.23) 

3,l,3',V 



t 



dd 



^dij'jj' 



77' 



E 

3,1,3', I 



(2.24) 

The expressions of the indirect d-d transfer integrals and 
the 3d-level energies thus obtained are equivalent to those 
obtained through the higher-order perturbational expan- 
sion with respect to the d-p and p-p transfers whose per- 
turbational processes contain not only one O 2p state but 
also more than one O 2p states. 

As a result, the "effective" multiband Hubbard Hamil- 
tonian derived from the multiband d-p model has the 
form; 



on— site? 



with 



(2.25) 



(2.26) 



Tj-mH 
^tdd 



= E 



<^^,,,^,4,d..y. + h.c.,(2.27) 



-ffon-site = Hu + H^' + Hj + Hj/ , (2.28) 

Here, -ff™^ expresses the energy levels of the 3d orbitals 
in the crystal field from the ligand oxygens. is the 

d-d supertransfer term. 

The advantages of the above formalism are as follows: 

• In the path-integral formalism, we can evaluate the 
amplitudes of the d-d transfers by treating the per- 
turbational processes which contain more than one 
O 2p states also as the intermediate states system- 
atically. 

• The obtained formulas of the d-d transfer and the 
crystal field terms are expressed by using parameters 

of the multiband d-p model. Since the values of these 
parameters have been well evaluated both theoreti- 
cally and experimentally, we can estimate the ampli- 
tudes of the d-d. transfer integrals and the magnitude 
of the crystal field quantitatively. 

By applying the second-order perturbational expan- 
sion to thus obtained multiband Hubbard Hamiltonian, 
wc derive an effective Hamiltonian on the subspace of 
states only with singly-occupied t2g orbitals at each Ti 
site in the insulating limit. We carry out the following 
expansion in terms of d-d transfers and j/U: 

1 



(2) 
eflF 



mH 

tdd 



Eq — {Hu+u' + Hj+j') 



HI 



mH 

tdd 



mH ^ 

tdd TP TT 
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j+j' 



j+3' 



Eq — Hu+u' 



-H. 



where 



Hu+u' — Hu + Hu' , 



j+3' 

(2.29) 

(2.30) 
(2.31) 

When one of the threefold t2g orbitals is occupied at 
each site, we can describe the electronic state using a 
spin-1 operator, which we call the pseudospin r. We 
express the occupied xy orbital by a quantum number 
Tz = —1, the yz orbital by = and the zx orbital 
by = +1. In this subspace, we rewrite the 2>d electron 
operator d^im'a'dima as Am,m'Ba^a' in terms of S and 
r depending on the indices m, m', a and cr'. Here, the 
indices m and m' run over the xy yz and zx orbitals. 
The symbol a and a' denote the spin of the 3d electron. 
The definitions of Am,m' and -Bo-.u' are as follows, 



1/2 t,(t,-1) 




for 


(m. 


to' 


= (1,1) 






for 


(to. 


to' 


= (2,2) 


1/2 T,(r^ + 1) 
1/V2 T+/3-^ 




for 


(to. 


to'" 


= (3,3) 




for 


(m. 


to' 


= (1,2) 


1/2 T+T+ 




for 


(to. 


to') 


= (1,3) 






for 


(to. 


to' 


= (2,1) 


1/^2 T+f3° 




for 


(m. 


to' 


= (2,3) 


1/2 T-T" 




for 


(m, to/ 


= (3,1) 


l/y/2 T-(3+^ 




for 


(to. 


to' 


= (3,2) 



(2.32) 



where 



and 



1/2 +5^ for 


(^, 


^') = 


(T,T) 


1/2 -Sz for 


(^, 


<^') = 


(i,i) 


S+ for 


(^, 


a') = 


(T,i) 


S- for 


(^, 


a') = 


a,T), 


T+ = Tj; + iTj;, 


r~ 







<S'~'~ iSa; + iSy, S Sx i-Sy 



(2.33) 



(2.34) 



(2.35) 



Here, the numbers to (to')=1, 2, and 3 correspond to the 
xy, yz and zx orbitals, respectively. 

We express the 3d electron operators in terms of S 
and r in accordance with the above rules to arrive at 
the effective spin-pseudospin Hamiltonian; 

frmH I Tj 



(2.36) 



The first term H™^ denotes the crystal field from the O 
ions acting on the Ti t2g orbitals. Without any distor- 
tion of the TiOe octahedron, the crystal field (i?™^) has 



a cubic symmetry so that it does not lift the t2g degen- 
eracy. In the following calculations, in order to study the 
effects of the R crystal field, we examine the case of no 
octahedral distortion. Therefore, we eliminate this term 
since it gives only a constant contribution. The effects 
of the O crystal field due to the actual TiOe distortion 
will be examined later. The second term Hgr is the spin- 
pseudospin term. We introduce an additional term for 
the crystal field from the R cations {Hay.). Finally, we 
obtain the Hamiltonian: 



(2.37) 



It should be mentioned that our approach is appropri- 
ate for the systematic study on properties and mecha- 
nism of magnetic and orbital phase transitions since the 
origin of the stabilization of magnetic and orbital ordered 
states are attributed to the second-order perturbational 
energy gains with respect to the indirect d-d transfers. In 
addition, their phase transitions are caused by the com- 
petition of these energy gains. The energy gains are easily 
evaluated in our approach by estimating the anisotropic 
electron-transfer amplitudes and the energy levels of the 
3d orbitals. 

3. Results on LaTiOa 

In this section, we study effects of the crystal field 
due to the displacements of the La ions induced by the 
GdFeOa-type distortion in LaTiOs. With the GdPeOs- 
type distortion, the La ions shift mainly along the 
(l,l,0)-axis (namely, the 6-axis), and shghtly along the 
(1, — 1, 0)-axis (the a-axis).^^'^^'^^^ There are two kinds 
of LaO-planes (plane 1 and plane 2) stacking alternat- 
ingly along the c-axis ((0, 0, l)-axis) as shown in Fig. 4. In 
plane 1, the La ions shift in the negative direction along 
the 6-axis while they shift in the positive direction in 
plane 2. Consequently, the crystal field from the La ions 
is distorted from a cubic symmetry. For sites 1 and 2, the 
distances between the Ti ion and the La ions located in 
the ±(1, 1, l)-directions decrease while those along the 
±(1, 1, — l)-directions increase (see Fig. 4 (a)). On the 
other hand, the distances along ±(1, 1, — l)-directions de- 
crease while those along ±(1, 1, l)-directions increase for 
sites 3 and 4. (In Fig. 4 (b), the stacking of each site is 
presented.) The changes of the other Ti-La distances are 
rather small. Since the nominal valence of the rare-earth 
ion is 3+, we expect that the Ti 3d orbitals directed along 
the shorter Ti-La bonds are lowered in energy because 
of the attractive Coulomb interaction. Without any dis- 
tortions of the TiOg octahedra, the crystal field from the 
ligand oxygens has a cubic symmetry. 

At this stage, when we introduce the crystal field from 
the La cations {Hm) as a perturbation, the threefold 
degeneracy of the cubic-t2g levels splits. Here, for sim- 
plicity, we take Hri by assuming that there exist point 
charges with +3 valence on each La cation. The Coiilomb 
interaction between an electron on a Ti 3d orbital and 
surrounding La^+ ions is given by using a dielectric con- 
stant CTiLa as follows, 



v{r) 



^ CTiLal-Ri - 



(3.38) 
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Fig. 4. (a) Displacements of the La cations around a TiOe octa- 
hedron (site 1) in the (JdFeOs structure are presented by arrows. 
In LaO-plane 1, they shift in the negative direction along the b 
axis while in LaO-plane 2, they shift in the positive direction. 
The distances between the Ti ion and two La ions (gray circles) 
located in the ±(1, 1, l)-directions are decreased. As a result, the 
La ions generate a crystal field similar to the D-^j crystal field 
with [1,1,1] trigonal axis (see text). Note that rotations of the 
octahedra due to the GdFeOs-type distortion are not presented, 
(b) Stacking of the TiOe octahedra as well as shifts of the La 
ions are presented. Note that Ti and O ions are not explicitly 
presented in this figure. 



where Ri expresses the coordinate of the i-th La ion, 
and Zfi{=+3) is the nominal valence of the La^+ ion. 
We calculate the following matrix elements, 

{m\v{r)\m') = J dr^* ^^^^v{r)ip3d,m' (3.39) 



with 



81V30 



Zm 
ao 



3/2 



(3.40) 
(3.41) 



Here, the integer m(=— 2, 2) is the magnetic quantum 
number, o,o tlic Bohr radius, p = ^f^r, and .^m(=+3) 
is the nominal valence of the Ti^+ ion. The coordinates of 
the La ions are calculated by using the positional param- 
eters and the cell constants obtained by the x-ray diffrac- 
tion study. After transforming the basis from m to the 
t2g representations, we obtain the Hamiltonian Hm with 
xy-, yz- and zx-basis. These representations are defined 
in terms of the x-, y- and ^;-axes attached to each TiOe 
octahedron. On the basis of the analyses of thus obtained 
iJiji, we show in the following that the consideration of 
the La crystal field well explains the emergence of the 



t2g 







cubic R crystal field 

Fig. 5. Energy-level structure in the crystal field due to the La 
cations. The threefold degenerate cubic-t2g levels split into three 
nondegenerate levels. 



AFM(G) ordering and the puzzling experimental find- 
ings in LaTiOs. 

By diagonalizing Hjn, we can evaluate the 3d-level 
energies and their representations. The threefold degen- 
erate t2g orbitals split into three isolated levels as shown 
in Fig. 5. The energy levels Ai and A2 are 0.77/eTiLa eV 
and 1.58/eTiLa eV, respectively. Here, it should be noted 
that since it is difficult to evaluate the local screening ef- 
fect in a solid, it is hard to determine the value of effective 
dielectric constant exiLa without ambiguity. Therefore, in 
our calculation, cxiLa is left as a variable. 

The representations of the lowest levels at each site are 
specified by the linear combinations of the xy, yz and zx 
orbitals as follows, 

site 1; a\xy) + c\yz) + h\zx), 

site 2; a\xy) + h\yz) + c\zx), 

site 3; a\xy) — c\yz) — b\zx), 

site 4; a\xy) - h\yz) - c\zx). (3.42) 

where a^+6^+c^=l. Here, we note that there exists 
a mirror plane vertical to the c-axis, and the orbital 
structure has the same symmetry as the orthorhombic 
GdFeOa-typc distortion. 

For LaTiOa, the coefficients take a=0.60, 6=0.39 and 
c=0.69. The orbital wavefunctions at sites 1 and 2 are 
similar to '^{^V + yz + zx), which is the lowest orbital 
for the Dsd crystal field with [1,1,1] trigonal axis. On the 
other hand, the orbital wavefunctions at sites 2 and 4 are 
similar to -^{xy — yz — zx), which is the lowest orbital 
for the D^d crystal field with [1,1,-1] trigonal axis. Since 
there are four trigonal axes in an octahedron, there are 
several configurations of the trigonal axes at sites 1, 2, 3 
and 4. The orbital structure realized in the crystal field of 
Hfii is similar to that with the lowest orbitals in the D^d 
crystal fields with [1,1,1], [1,1,1], [1,1,-1] and [1,1,-1] 
trigonal axes at sites 1, 2, 3 and 4, respectively (see 
Fig. 6). In the previous study (Ref. 40), this trigonal-axes 
configuration was referred to as config. 3, and it was pro- 
posed that the orbital ordering in the D^d crystal fields 
with config. 3 strongly stabilizes the AFM(G) ordering, 
and well explains the puzzling experimental properties 
of LaTiO;j. In addition, recent dynamical mean- field cal- 
culations combined with first principles local density ap- 




Fig. 6. Orbital structure in the D^^i crystal fields with [1,1,1], 
[1,1,1], [1,1,-1] and [1,1,-1] trigonal axes at sites 1, 2, 3 and 4, 
respectively. The orbital wavefunctions at each site are {xy + 

yz + zx), -^{xy + yz + zx), -^{xy - yz - zx) and -^{^y - 
yz — zx), respectively. Note that tiltings due to the GdFeOa- 
type distortion are not presented explicitly. 

proximation (LDA+DJVIFT) show similar orbital orders 
with nondegenerate orbital ground state for LaTiOa and 
are consistent with the results obtained here.^^'^°) 

Here, we note that this orbital ordering likely induces 
the I?3d distortions of the unit TiOg octahedra with the 
trigonal axes of config. 3. Recently, such distortions were 
actually observed by diffraction measurements.^^'' One 
might argue that this TiOg distortion could occur by the 
usual Jahn- Teller mechanism. However, this is not true 
because the orbital degeneracy is already lifted by the R 
crystal field and the Jahn- Teller mechanism is not effec- 
tive if the GdFeOs-type distortion (tilting of the TiOg 
octahedra) is present. The observed TiOg distortion is 
simply induced by the orbital ordering due to the R crys- 
tal field. In this sense, the t2g degeneracy is lifted mainly 
by the R crystal field, and the TiOg distortion occurs 
subsequently. Our approach in which the R crystal field 
is taken first as the primary driving mechanism is based 
on the fact that the GdFeOa-type distortion is stabilized 
from much more high-energy origin of the R-0 covalency 
as is detailed in [2] of Sec. 5. 

We next discuss the stability of the magnetic state in 
the crystal field of Hpn. For this purpose, we employ the 
effective spin-pseudospin Hamiltonian introduced in the 
previous section. Substituting Hpn for the crystal field 
term i?cry., we calculate energies for several magnetic 
structures by applying the mean-field approximation. In 
Fig. 7, we plot the calculated energies as functions of the 
magnitude of Ai (=0.77/eTiLa)- The magnitude of Ai 
is varied because of the uncertainty due to exiLa- In the 
region of Ai >0.03 eV, the AF]VI(G) state is strongly 
stabilized relative to the other structures. 

In this region, the value of Ai is much larger than 
/cb2n so that the orbital occupation is restricted to the 
lowest orbitals irrespective of the spin structure. Now, 
we estimate the magnitude of the spin-exchange inter- 
action for each Ti-Ti bond in the subspace of singly- 
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Ai (eV) 

Fig. 7. Energies for several magnetic structures under the crystal 
field of Hfii are plotted as functions of Ai. 

occupied lowest orbitals. The spin-exchange constant J 
for a Ti-Ti bond is represented as J = — £^||)/2S'^ 
with and E-^^ being the mean- field energy gains for 
the Ti-Ti bond of tt" and fj-pairs, respectively. For 
LaTiOa, the values of J along the x-, y- and z-axes take 
as J3; = I8.5 meV, Jy=I8.5 meV and Jz=19.7 meV, re- 
spectively, which well reproduce the spin-wave spectrum 
with isotropic exchange constant of ~ 15. 5 meV. 

We have confirmed that the values of J as well as the 
orbital wavefunctions obtained by the present calculation 
hardly change (within the accuracy of a few %) even if 
the next nearest neighbor O and R ions are taken into 
account, which suggests that the present result should 
be well established even with the ions further than the 
nearest neighbors since in reality, the ion-ion interaction 
is strongly screened as the distance increases. The ro- 
bustness of the above results will be discussed in more 
detail in Sec. 5. 

Here, we further note that the crystal field from the 
La cations has two origins. One is the Coulomb potential 
from charged La ions as we have studied above, and the 
other is the hybridization between the Ti 3c? orbitals and 
unoccupied orbitals on the La ions. We next examine the 
latter effects for LaTiOs. We construct the Hamiltonian 
for hybridization between the Ti 3c? and La 5c? orbitals 
{Hr2) by using the second-order perturbational expan- 
sion in terms of the transfers between the Ti 3d and La 
bd orbitals [t'^'^). The expression of the matrix element 
of Hr2 is 

j-dd j-dd 

{m\Hn2\m') = -Y,^^^^Y^- (^.43) 

2,7 

Here, the indices m and m' run over the cubic-i2g rep- 
resentations, xy, yz and zx. The symbols i and 7 are 
indices for the eight nearest-neighbor La ions and the 
fivefold La bd orbitals, respectively. The symbol As^^ de- 
notes the characteristic energy difference between the 
Ti t2g and La hd orbitals. The transfer integral i^"* is 
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Table I. Energy-level structures and representations of the lowest 
orbitals for Hjn and Hji2- 







Hri 


Hr2 


Ai 


(eV) : 


0.77/eTiLa 


0.08 


Aa 


(eV) : 


1.58/eTiLa 


0.14 


a : 




0.60 


0.62 


b : 




0.39 


0.44 


c : 




0.69 


0.65 



given in terms of Slater-Koster parameters Vdda, VddTr 
and VddS- It is assumed that these parameters are pro- 
portional to d^^ with d being the Ti-La bond length. ^'^^ 
On the basis of the analyses of LDA band structure^"'^^ 
and the photocmission spectra, we fix these parameters 
as Vdda = -1.04 eV, Vdd-. = 0.56 eV and VddS = eV for 
the Ti-La bond length of 3.5 A, and A^d = 3.6 eV. 

By diagonalizing thus obtained i?H2j we have evalu- 
ated the energy levels and their representations (see Ta- 
ble I). In the obtained energy level structure, the three- 
fold t2g levels split into three nondegenerate levels, which 
is similar to that of Hm . Since the representation of the 
lowest orbital of Hr2 is also similar to that of Hjn, the 
level sphtting of H]ii+Hji2 is well expressed by the sum 
of these two contributions. Considering the fact that 
AFM(G) state is stabilized in the region of Ai >0.03 eV, 
and Ai for Hji2 is ~0.085 eV, the crystal field due to the 
hybridization between the Ti 3d and La 5d orbitals {Hr2) 
alone turns out to stabilize the AFM(G) spin structure 
strongly. 

In addition, the value of Ai is sufficiently large as com- 
pared with the coupling constant of the LS interaction 
in Ti3+ ((^=0.018 eV).^^) Consequently, the crystal field 
from the La cations dominates over the LS interaction, 
resulting in the quenched orbital moment. Now, the ex- 
perimentally observed reduction of the magnetic moment 
is attributed to the strong itinerant fluctuation near the 
metal-insulator phase boundary instead of the LS inter- 
action as discussed in the literature.*"-' We will discuss 
this issue in § 5 in more detail. 

4. Results on HTiOg (J?=La, Pr, Nd and Sm) 

In this section, we study the effects of the crystal field 
from the R ions in /JTiOs as functions of the GdFeOs- 
type distortion (in other words, as functions of the size of 
R ion). When the TiOe octahedra forming the perovskite 
lattice are distorted, the crystal field from the ligand 
oxygens is considered to be responsible for the electronic 
structure. However, with small or no atomic distortion of 
the octahedra, we expect that the crystal field from the 
R ions sufficiently lifts the t2g and Cg degeneracy of the 
3d state, and plays a substantially important role on the 
stability of the magnetic and orbital states. The magni- 
tude of the local distortion of the TiOe octahedron is ex- 
pressed by the standard deviation from the mean for the 
three inequivalent Ti-0 bond length. In Fig. 8, we plot 
the standard deviation (tr) obtained by the x-ray diffrac- 
tion study^^-* as a function of the Ti-O-Ti bond angle. 
The value for PrTiOa is obtained by the interpolation of 
LaTiOa and NdTiOa. Focusing on the compounds with 
relatively small a, we have studied relations between the 
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Fig. 8. Standard deviation from the mean for three inequivalent 
Ti-O bond length is plotted as a function of the Ti-O-Ti bond 
angle. 

GdFeOa-type distortion and the R crystal field system- 
atically. We again derive the crystal field Hamiltonians 
based on the point charge model (i?i?,i) for RT'iO^ with 
R being Pr, Nd and Sm by using the coordination pa- 
rameters, and calculate the wavefunctions of the lowest 
orbitals. Here, the value of a for SmTiOa is rather large 
as compared with the other three compounds of i?=La, 
Pr and Nd, which indicates relatively strong distortion 
of the TiOe octahedra. This means that the crystal field 
from the O ions in SmTiOa is deviated from the cu- 
bic symmetry, which would strongly affect the electronic 
structure in this compound. Therefore, we have to take 
effects of both Sm and O crystal fields into account. Thus, 
we first study the effects of the Sm crystal field {Hm). 
Then, the effects of the distortion of the TiOe octahedra 
and the O crystal field will also be studied later. 

First, we study the orbital structures in the R crystal 
fields. In LaTiOs with the smallest GdFeOs-type dis- 
tortion, the closest two R ions are located nearly along 
the trigonal directions, and the orbital wavefunctions 
are approximately expressed by -^[xy + yz + zx) and 

-^(a;?/ — yz — zx) as discussed in the previous section. 
The weights of the xy, yz and zx orbitals are almost 
the same in this compound. As the GdFeOa-type dis- 
tortion increases, the orbital wavefunctions should differ 
from those of LaTiOa since the locations of the R ions 
are deviated from the trigonal directions. We expect that 
with increasing GdFeOs-type distortion, the component 
of the zx orbital gradually decreases at sites 1 and 3, 
while the yz-component decreases at sites 2 and 4. This 
expectation is based on the following discussion. Tilting 
of the TiOe octahedra in the GdFeOa-type distortion is 
approximately expressed by rotations around the h- and 
c-axes (see also Fig. 9). The rotation around the c-axis is 
clockwise at sites 1 and 3, and counterclockwise at sites 
2 and 4. Therefore, at sites 1 and 3, distances between 
the zx-plane and the closest two R cations increase with 
increasing distortion, while those for yz-plane increase at 
sites 2 and 4. As a result, the attractive potential between 
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Fig. 9. Rotations of the TiOg octahedra around the c-axis in the 
GdFeOa structure are presented for sites 1-4. The closest two _R 
cations are also presented by gray and white large circles. 




148 150 152 154 156 158 
Ti-O-Ti bond angle (deg) 



Fig. 10. Orbital wavefunctions at each Ti site are presented as 
functions of the Ti-O-Ti bond angle. The dashed line indicates 
1/V3. When the bond angle decreases, the orbitals continuously 
approach the character of YTiOa given by Eq. (4.44). 



the R cations and the zx orbitals are reduced at sites 1 
and 3, while those for the yz orbitals are reduced at sites 
2 and 4, resulting in the decrease of the corresponding 
orbital occupations. 

In this discussion, only the rotations around the c-axis 
are taken into account. However, in the actual situation. 



the rotations around the 6-axis also exist. Moreover, the 
locations of the R ions are different from the ideal po- 
sitions assumed in the above discussion. Thus, for more 
precise estimate, it is necessary to study the orbital struc- 
tures by constructing the crystal field Hamiltonian Hm 
in which the experimentally measured position parame- 
ters of the R and O ions are taken into account. 

The lowest orbitals of Hj^i for these compounds are 
also expressed by Eq. (3.42). In Fig. 10, we show the co- 
efficients a, b and c for the lowest orbitals of Hjf i at each 
Ti site as functions of the Ti-O-Ti bond angle. With rel- 
atively small GdFeOa-type distortion, the orbital wave- 
function has a form similar to that of the lowest orbital 
in the D^d crystal field, namely '^{^V + y-^ + zx) or 
-^{xy — yz — zx). As the GdFeOa-type distortion in- 
creases, the component of the zx orbital actually de- 
creases at sites 1 and 3, while the yz-component de- 
creases at sites 2 and 4, which is consistent with the 
above discussion. 

The reductions of these orbital occupations indicate 
that with increasing GdFeOa-type distortion, the crystal 
field from the R cations tends to progressively stabilize 
the orbital structure in which sites 1, 2, 3 and 4 are 
occupied by the orbitals as follows: 



site 



site 



site 



site 



1; 



\yz)), 



1 

vl 

1 

V2 



{\xy) + \zx)), 



{\xy) - \yz)), 



v)). 



(4.44) 



In this orbital structure, there exists a mirror plane 
which is vertical to the c-axis. In fact, this orbital struc- 
ture is the same as the orbital ordering which has been 
observed in YTiOa with large d-type Jahn- Teller distor- 
tion by several experiments such as NMR (Ref. 36,37), 
polarized neutron scattering (Ref. 64,65), and resonant 
x-ray scattering (Ref. 66). Further, a recent resonant 
x-ray scattering study shows that the orbital state in 
SmTiOs has a twofold symmetry similarly to YTiOs and 
GdTiOs, which is consistent with the present result. 
It is interesting to note that, in our result, the character 
of the orbital symmetry in YTiOs is retained but quan- 
titatively and continuously decreases toward the type of 
LaTiOs. The detection of the orbital ordering in each 
compound predicted here can determine the validity of 
our theory. 

Now let us discuss the relations between the orbital 
state here obtained and the experimentally observed 
Jahn- Teller distortion in YTiOa and GdTiOa . In the pre- 
vious section, we have mentioned that in LaTiOa, the 
orbital ordering stabilized due to the La crystal field fa- 
vors the stacking of the trigonally distorted TiOe octahe- 
dra with [1,1,1] and [1,1,-1] trigonal axes alternatingly 
along the c-axis, which was actually detected experimen- 
tally. On the other hand, YTiOa and GdTiOa with strong 
GdFeOa-type distortion show the d-type Jahn- Teller dis- 
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Fig. 11. d-typc Jahn- Teller distortion. The threefold degenerate 
t2g levels split into twofold degenerate lower levels and a non- 
degenerate higher level at each Ti site. The ways of splitting at 
each site are also shown. 



tortion in which the elongated axes of the octahedra are 
paraUel along the c-axis as shown in Fig. 11. As a result 
of this distortion, the xy and yz orbitals are lowered in 
energy at sites 1 and 3, while at sites 2 and 4, the xy and 
zx orbitals are lowered. In fact, the orbital wavefunctions 
in Eq. (4.44) are represented by linear combinations of 
the twofold lower orbitals in the d-type Jahn- Teller dis- 
tortion. This suggests that the d-type Jahn- Teller dis- 
tortion realized in i?Ti03 [R —Y and Gd) is favored by 
the orbital ordering stabilized by the R crystal field. In 
fact, a recent multiband d-p model study has shown that 
with strong GdFeOa-type distortion, the covalency be- 
tween the O and R ions induces the d-type Jahn- Teller 
distortion.^*) We stress that these two effects may co- 
operatively work to stabilize this distortion. It is inter- 
esting to point out that in the titanates, the R ions in 
the GdFeOa-type structure may control the octahedral 
distortions. 

Here, we note that the orbital structure in SmTiOs 
may be slightly deviated from the orbitals presented here 
because of an additional effect of the O crystal field with 
relatively large distortion of the TiOe octahedra. More 
concretely, we expect that in SmTiOa, the occupation of 
the xy orbital is somewhat increased at all Ti sites as 
will be discussed later. 

We also calculate the spin-exchange constant J by as- 
suming that the orbital occupation is restricted to the 
lowest orbital at every Ti site. In Fig. 12, we plot the 
values of J^, Jy and Jz as functions of the Ti-O-Ti bond 
angle. The AFM(G) spin structure is also reproduced in 
i?Ti03 with R being Pr, Nd and Sm. In LaTiOs with 
the smallest GdFeOa-type distortion, a nearly isotropic 
spin coupling is realized. As the GdFeOa-type distortion 
increases, the spin exchange gradually decreases. This is 
consistent with the gradual decrease of Tn as R goes 
from La, Pr, Nd to Sm. In this calculation, the decrease 
of Jx and Jy is rather steep relative to that of Jz , result- 
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Ti-O-Ti bond angle (see text). 



ing in the slightly anisotropic spin couplings. However, 
these subtle anisotropies can possibly be diminished or 
affected when we introduce the crystal field from the O 
ions as will be discussed later. It should be noted that 
these anisotropies will not necessarily be observed exper- 
imentally. 

Here, we note that Tn for SmTiOa is about 50 K, 
which is decreased to 40% of Tn for LaTiOa. However, 
the calculated spin-exchange constant is not so strongly 
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depressed, and the value of Jx for SmTiOa is about 70% 
of that for LaTiOs. This is puzzHng when we consider the 
fact that the model calculation based on the 3D Heisen- 
berg model on the cubic lattice show that Tn is propor- 
tional to J.s"^) 

In order to solve this problem, we study effects of the 
distortion of the TiOg octahedra. In the above analyses, 
we have examined the effects of the GdFeOa-type distor- 
tion by assuming the undistorted TiOg octahedra as well 
as the cubic O crystal fields since the generic GdFeOs- 
type distortion induces only tiltings and rotations of the 
octahedra, and does not distort them. However, in ad- 
dition to the GdFeOa-type distortion, distortions of the 
octahedra themselves are relatively strong in SmTiOa so 
that it is necessary to examine the effects of the O crys- 
tal field deviated from the cubic symmetry in order to 
clarify the electronic structure in this compound. 

We derive the crystal field Hamiltonians for the repul- 
sive Coulomb potential between an electron on the Ti 
orbitals and the ligand ions (i?Li) by assuming the 
form of the Coulomb interaction as 



v{r) 



2e' 



exiol-Ri 



(4.45) 



where exio is a dielectric constant, and Ri denotes the 
coordinates of the i-th O ions. 

In addition, the crystal field Hamiltonians for hy- 
bridization between the Ti and O 2p orbitals {H]^2) 
is also derived. The expression of the matrix element of 
Hl2 is 



rpa ,pa 

I 17-7- I /\ X ^ ra\il m' :il 

{m\HL2\m) =2_^' 



A 



(4.46) 



Here, the indices m and m' run over the cubic-t2g repre- 
sentations, xy, yz and zx. The symbols i and I are indices 
for the six ligand O ions and the threefold O 2p orbitals, 
respectively. The symbol A stands for the charge-transfer 
energy. The transfer integral between the Ti 3d and O 
2p orbitals {tP''-) is given in terms of Slater-Koster pa- 
rameters Vpda, Vpd-K- These parameters are calculated by 
assuming that they are proportional to d~^'^ with d be- 
ing the Ti-0 bond length.^^^ 

In Fig. 13, we plot relative energies of each level for 
Hli (upper panel) and those for (lower panel). In 
the compounds with small GdFeOs-type distortion, the 
differences of the energies for each level are rather small 
whereas in largely distorted GdTiOs and YTiOs, one 
orbital is higher in energy relative to the other two or- 
bitals. In GdTiOa and YTiOa, the xy and yz orbitals 
are strongly lowered in energy at sites 1 and 3, and the 
xy and zx orbitals at sites 2 and 4 with the large d-type 
Jahn- Teller distortion. 

In the case of moderately distorted SmTiOa, the O 
crystal field lowers the xy orbital at every Ti site, and 
both yz and zx orbitals are higher in energy. As a result, 
the distortions of the TiOg octahedra in SmTiOs favor 
the xy-occupation at every site. In this compound, the 
crystal field from the Sm ions and that from the O ions 
are strongly competing, and the low-energy orbital-spin 
structure is determined by this competition. In SmTiOa, 
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Fig. 14. (a) Energy-level scheme of the Sm crystal field Hamilto- 
nian Hjn for SmTiOs. (b) Energy-level scheme of the O crys- 
tal field Hamiltonian Hi^i for SmTiOs. (c) The lowest-orbital 
wavcfunctions in the crystal fields from both Sm and O ions 
{Hri+Hli) as functions of Ao/Asm- (d) Spin-exchange con- 
stants along the x-, y- and z-a-KCS as functions of Aq/Ashj. 



the Sm crystal field [Hm) tends to lift the t2g degen- 
eracy into three isolated levels as shown in Fig. 14 (a). 
The magnitude of splitting Asm is 1.44/eTiSm eV, and 
the lowest-orbital wavcfunctions at each site are repre- 
sented by Eq. (3.42) with a =0.63, b = 0.25 and c = 0.74. 
On the other hand, if only the O crystal field {Hli) ex- 
ists, the t2g degeneracy splits into lower xy orbital and 
higher yz and zx orbitals with Aq of 0.20/eTio eV as 
shown in Fig. 14 (b). Because of the uncertainty due to 
CTiSm and exiOj it is difficult to estimate the strength of 
these two crystal fields qualitatively. Then, we study the 
orbital structure and the spin coupling by varying the 
ratio of the strength of these two crystal fields (in other 
words, by varying Ao/Asm)- The lowest-orbital wave- 
functions are also expressed by Eq. (3.42) even when 
both O and Sm crystal fields (i?_Ri and Hli) exist. In 
Fig. 14 (c), we plot the coefficients a, b and c as func- 
tions of Ao/Asm- This figure shows that the xy-orbital 
occupation increases with increasing Ao/Asm, namely, 
with increasing strength of the O crystal field. 
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The increased .Ty-occupation is expected to cause a 
two-dimensional anisotropy in the spin coupling in which 
a strong AFM coupling is realized in the a6-plane while 
along the c-axis, the AFM coupling is rather weak. In 
Fig. 14 (d), the spin-exchange constants J^, Jy and 
Jz are plotted as functions of Ao/Asm- In this figure, 
the value of Jz rapidly decreases while Jx and Jy in- 
crease as Ao/Ag„i increases. In particular, in the re- 
gion of Ao/Asni > 1, the value of becomes almost 
zero while and Jy take relatively large value, and 
the ratio Jz/Jx becomes smaller than 0.13. This means 
that when the strength of the O crystal field is compa- 
rable to, or is larger than the Sm crystal field, an al- 
most perfect two-dimensional spin coupling is realized. 
These facts indicate that the increased xy-occupation 
due to the O crystal field actually causes the two- 
dimensional anisotropy in the spin coupling. The ob- 
served depression of Tn in SmTiOa may be attributed 
to this two-dimensional anisotropy. We predict that this 
two-dimensional anisotropy may be observed in SmTiOs. 

In addition, the crystal field from the Sm ions and 
that from the O ions make opposite contributions to the 
anisotropy of the spin coupling. As shown in Fig. 14, 
without the O crystal field (namely, Ao/Asm=0), the 
Sm crystal field favors the anisotropy in which spin ex- 
change along the c-axis is slightly larger than that in the 
a6-plane. On the other hand, the O crystal field favors 
the in-plane anisotropy. Figure 14 shows that the two- 
dimensional anisotropy due to the O crystal field domi- 
nates over the small contribution to the anisotropy from 
the R crystal field even in the region of small Ao/Asm- 

In SmTiOs, we expect that both crystal field from the 
Sm ions and that from the O ions play important roles 
on the electronic state, and the orbital-spin structure in 
the Sm crystal field is strongly affected by the O crystal 
field because of the relatively strong TiOe distortion. It is 
essentially important to study the orbital state by quan- 
titative estimate of the crystal field from both O and Sm 
ions to elucidate the orbital-spin structure of SmTiOa, 
which is left for further study. 

In the compounds of R =Pr and Nd, if we neglect the 
effects of the octahedral distortion, the calculated spin- 
exchange constants show very small anisotropics in which 
the AFM spin coupling along the c-axis is slightly larger 
than that in the a6-plane as shown in Fig. 12. However, 
even when the magnitude of the level splitting is consid- 
erably small, the O crystal field in PrTiOa and that in 
NdTiOs also favor the ajy-occupation, which generates 
the AFM coupling in the a5-plane. Under these circum- 
stances, the slight anisotropy favored in the R crystal 
field should be diminished or weakened by the subsequent 
and opposite contribution from the O crystal field. Thus, 
we expect that the AFM coupling in PrTiOa and that in 
NdTiOs can be rather isotropic. The quantitative esti- 
mates of the competition between R and O crystal fields 
in these compounds are left for further studies. 

We finally mention that when the crystal field from 
the O ions due to the TiOg distortion is superimposed, 
the dominance of the R crystal field is taken over by the 
O crystal field at ii=Sm. Through severe competition 
at i?=Sm, the O crystal field deviated from the cubic 
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symmetry due to the d-type Jahn- Teller distortion comes 
to control the orbital-spin structure at i?=Gd and Y. 
We note that in the present analyses, we have focused 
on the effects of the R crystal fields on the cubic-t2g 
degeneracy, and the Jahn- Teller distortion is not taken 
into account so that the orbital-spin structure in the FM 
region with strong d-type Jahn-Teller distortion has not 
been reproduced. Therefore, a recent puzzling finding of 
isotropic ferromagnetism in YTiOs^*^ is not resolved in 
this study and is left for further study. 

5. Discussion 

In this section, we will discuss five points on our pro- 
posal. Before the detailed discussion, we first summarize 
these points in the following: 

[1]: Our theory succeeded for the first time in explaining 
puzzling experimental results in a unified way. The avail- 
able experimental results will be discussed in the light of 
our theory. 

[2]: A universal mechanism for controlling electronic 
structure of the perovskite-type transition-metal oxides, 
which has been overlooked is proposed. The mechanism 
here proposed is substantially different from the usual 

Jahn-Teller mechanism. 

[3]: We will also discuss the accuracy of the numerical 
results and the reliability of our predictions about orbital 
structures as well as the spin-coupling isotropy. 
[4]: Based on the present results together with the pre- 
vious ones,^^'^^'^^'^'^'^^^ the overall features of the mag- 
netic phase diagram for the titanates will be discussed. 
[5]: The observed magnetic moment of 0.45 /Ub in 
LaTiOs^^^ has been sometimes considered to be small 
by comparing with the value of ~0.85 /^b in the 
spin-wave theory on the 3D Heisenberg model. We will 
discuss the origin of this reduction. 

[1]: In the present paper, we have shown that a very 
simple mechanism of the R crystal field here, we 
proposed and never been considered seriously before, 
accounts for available experimental results with no 
contradiction. Since the perovskite titanates attract an 
enormous amoiint of interest, a number of experiments 
have been done so far. However, there has been no 
theory which can explain these experimental results in 
a unified way. As discussed in the present paper, the 
following experimental results are explained within our 
theory for the first time. 

(1) : Emergence of the AFM(G) state 

The emergence of AFM(G) state in LaTiOs in which 
the t2g orbitals have been considered to be degenerate, 
has long been puzzling. We have shown in this paper 
that the AFM(G) state is realized with quenched orbital 
degrees of freedom due to the crystal field from the La 
ions. Further, the analyses of the R crystal fields in 
iiTiOs with R =Pr, Nd and Sm have also reproduced 
the AFM(G) ordering. 

(2) : Neutron scattering 

Recent neutron scattering experiment showed the spin- 
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wave spectrum characterized by isotropic spin coupling 
with J of 15.5 meV as well as a small spin gap for 
LaTiOs.^^) These characteristics are well reproduced 
in our picture. Since the lowest orbital state in the La 
crystal field is approximately represented by the linear 
combination of the xy, yz and zx orbitals with the same 
weights, nearly isotropic; elec;tron transfers along x, y 
and z directions are realized. The isotropic spin coupling 
is caused by this isotropy of transfer amplitudes. In 
addition, we have shown that the crystal field splitting 
is much larger than the coupling constant of the LS 
interaction, and the resultant quenched orbital moment 
leads to the observed small spin gap. 

(3) : Behavior of Tn in the magnetic phase diagram 
According to the experimentally obtained magnetic 
phase diagram,^^'^^'^^'^^^ Tn monotonically decreases 
as R goes from La, Pr, Nd to Sm with increasing 
GdFeOs-type distortion. The calculation of the spin 
coupling constant J based on our theory has well 
reproduced this behavior. 

(4) : Resonant x-ray scattering 

Recent resonant x-ray scattering experiment showed 
the evidence of an orbital ordering in the AFM(G) 
compounds of iJTiOs (R=La, Pr, Nd and Sm), and 
that the symmetry of the orbital ordering is the same as 
that in the FM compounds YTiOa and GdTiOs.^^) Our 
calculation of the orbital states based on our theory well 
accounts for this result. 

(5) : X-ray diffraction measurement 

We have shown that in LaTiOs, the crystal fields from 
the La cations stabilize an orbital ordering in which 
-^{xy + yz + zx), -^{xy + yz + zx), -^^(xy - yz - zx) 

and -^{xy — yz — zx) are approximately occupied at 
four Ti sites in the orthorhombic unit cell in an alter- 
nating fashion along the c-axis. This orbital ordering 
is expected to induce nearly trigonal distortions of the 
TiOe octahedra with trigonal axes of [1,1,1], [1,1,1], 
[1,1,-1] and [1,1,-1] at each Ti site. These distortions 
were actually observed by recent x-ray diffraction mea- 
surements.^*^ In addition, we have pointed out that the 
(i-typc Jahn- Teller distortions in YTiOs and GdTiOs is 
also favored by the orbital orderings expected in the R 
crystal field. 

(6) : NMR 

A NMR spectrum for LaTiOs proved to be well repro- 
duced if we take the orbital ordering model predicted in 
our theory. 

We note that there exists no other theory which 
explains these experimental results in such consistent 
ways, and that our theory succeeded in the explanation 
for the first time. 

[2]: We also emphasize the importance of our proposal 

as a universal mechanism determining the electronic 
structure in the perovskite-type transition-metal oxides. 



Through the present study, we have pointed out for the 
first time that in the perovskite-type transition-metal 
oxides, a crystal field from the R ions caused by the 
GdFeOa-type structure substantially lifts the Sd orbital 
degeneracy. This mechanism of lifting the orbital degen- 
eracy competes with the usual Jahn- Teller mechanism, 
and eventually plays an important role on controlling 
the orbital-spin structures. 

Our proposal includes two important points. One point 
is the significance of the R ions. The perovskite-type 
transition-metal oxides have often been studied by em- 
ploying some models which include effects of the strong 
electron correlations such as the Hubbard model, d-p 
model, t-J model and so on by assuming only the crystal 
fields on 3d orbitals from the surrounding octahedron. 
In these models, effects of the crystal field due to the 
R ions have been usually overlooked. However, we have 
shown that it is important to consider the roles of the 
R ions when we study the electronic structure of these 
compounds. 

Another important point is a new role of the GdFeOs- 
type distortion. In the GdFeOs-type distortion, TiOg oc- 
tahedra tilt alternatingly, which decreases the M-O-M 
bond angle from 180°. This distortion has been focused 
from the aspect of the bandwidth control since the de- 
crease of the M-O-M angle significantly reduces transfer 
integrals between neighboring transition-metal ions me- 
diated by the oxygen ions. On the other hand, an octa- 
hedron itself does not atomically distort in the generic 
GdFeOa-type distortion so that the O crystal field main- 
tains a cubic symmetry. Therefore, it had been consid- 
ered that the distortion does not lift the t2g and Cg degen- 
eracies. However, our result indicates that in reality, the 
distortion necessarily lifts the degeneracy through gener- 
ating the R crystal field, and determines the orbital-spin 
state. This direct effect on the electronic structure has 
been overlooked so far. 

Here, it should be noted that this mechanism is not 
a usual Jahn- Teller one. In the Jahn- Teller systems, 3d 
electrons on the transition-metal ions and the lattice 
degree of freedom strongly couple to each other, and the 
lattice spontaneously distorts to lift the 3d degeneracy 
and lower the electronic energy. On the other hand, 
the GdFeOs-type distortion is caused by interactions 
between the R ions and O ions. More concretely. Wood- 
ward showed that the octahedral tilting is induced by an 
energy gain in R-0 covalency bonding in the perovskite 
compounds with i?-site mismatching, and subsequent 
i?-site shifts occur by R-0 ionic interactions.^^) In 
Ref. 69, it has been shown that the displacement of the 
R ions is a direct consequence of the GdFeOs-type tilting 
through R-0 coupling. We note that this distortion 
exists irrespective of the 3d electronic state, and has 
nothing to do with any spontaneous lift of the electronic 
degeneracy in contrast with the Jahn- Teller mechanism. 
The lifting of the degeneracy of the Ti 3d orbitals occurs 
just as a consequence of the crystal field from the R 
ions. In this sense, the mechanism we have proposed 
is not a Jahn- Teller but an overlooked one, which is 
inherent and universal with GdFeOs-type distortion. 



16 J. Phys. Soc. Jpn. 

[3]: In the present study, we have examined the point 
charge model with the nearest neighbor 8 R ions in the 
case of degenerate t2g orbitals in the cubic crystal field 
from the nearest neighbor 6 O ions. The calculation of 
the spin-exchange constant well reproduced the experi- 
mentally measured isotropic superexchange parameters 
in LaTiOs. Now, let us discuss the adequacy of the ap- 
proximation used and the robustness of the predicted 
isotropy. 

First, we would like to stress that the present results 
would not be changed even if further ions would be taken 
into account. Our approximate treatment with the near- 
est neighbor ions actually provides correct results not 
only qualitatively but also quantitatively. 

In order to confirm this, we have calculated the 
t2g splittings, the orbital wavefunctions and the spin- 
exchange constants by using a point charge model with 
additional next-nearest-neighbor 24 R and 24 O ions. 
The differences between thus obtained values and those 
obtained by the point charge model with only nearest 
neighbor ions are considerably small (within a few %) 
as shown in the following: 

With only nearest neighbor O and R ions: 

• t2g splitting (Ai): 0.77/eLaTi 

• wavefunction: 0.605a;y -|- 0.390y^; -|- 0.695^;a; 

• spin exchange: Jab =18.5 meV, Jc =19.7 meV, 
With additional next nearest neighbor O and R ions: 

• t2g splitting (Ai): 0.75/eLaTi 

• wavefunction: 0.608a;y -|- 0.396yz + 0.688^a; 

• spin exchange: Jab =18.9 meV, Jc =19.4 meV. 

These results are calculated by assuming the same di- 
electric constant exiLa for both Coulomb repulsion from 
the nearest neighbor ions and that from the next nearest 
neighbors. In this sense, strong screening effects in a solid 
are neglected so that the differences are overestimated. 
Since the screening effects strongly depress the Coulomb 
interaction as the ion- ion distance increases, we expect 
that the difference should be much smaller if we take the 
screening properly into account. 

Here, we note that the purpose of the present paper 
is not a quantitatively precise estimate of the potential 
from the surrounding ions, but the proposal of a qualita- 
tive and transparent picture for the novel mechanism in 
which R cations lift the t2g degeneracy and determine the 
orbital-spin structure in the perovskite titanates. Never- 
theless, it should be emphasized that the point charge 
model with the nearest neighbor R and O ions here we 
employed gives sufficiently correct results also quantita- 
tively. 

In addition, in the present paper we have shown that 
the hybridization between the La 5d and Ti 3d orbitals 
works cooperatively with the Coulomb potential from 
the nearest neighbor R ions on the stabilization of the 
above orbital ordering and the AFM(G) state, and the 
hybridization alone has proved to strongly stabilize the 
orbital-spin structure. This fact also supports that the 
consideration of the Coulomb potential from further R 
ions hardly changes the results. 
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We also note that the predicted isotropy is robust 
even when we consider the slight TiOe distortion, which 
is actually realized in LaTiOa. In our analyses, the TiOe 
octahedron is assumed to be undistorted. We consider 
that the actual orbital state may be slightly deviated 
from that calculated here because of the small TiOe 
distortion. Actually, in the paper by Cwik et al., by 
considering the small octahedral distortion observed 
in their experiment, the orbital state is calculated as 
Q.n{-^{zx+yz))+Qm&xy = QMQxy+Q.b^yz + Q.bAzx. 

This is also similar to the -^[xy + yz + zx) orbital, 
but is slightly different from our result. However, in the 
previous study, we have evaluated the spin-exchange 
constant when the perfect -^{xy + yz + zx) orbital is 
occupied, and we have shown that the spin coupling is 
also almost isotropic and the anisotropy {Jab/ Jc) is less 
than 10 %. This result together with the present one 
suggests that the isotropy is hardly affected by the slight 
deviation of the orbital state from -^{xy + yz + zx). 
Therefore, we conclude that the isotropy predicted in 
our study is robust. 

[4]: We next discuss the experimentally obtained mag- 
netic phase diagram^^' based on the present re- 
sults together with the previous ones^'^' (see 
also Fig. 15). In LaTiOs with the smallest GdFeOa-type 
distortion, since the atomic distortion of the TiOe oc- 
tahedra is rather small, the crystal field from the La 
ions dominantly determines the Ti 3d state. As a result, 
LaTiOs exhibits a AFM(G) ordering with the lowest or- 
bital occupation in the nearly-trigonal La crystal field. 
The AFM(G) ordering due to the R crystal field is also 
realized in the compounds of Pr, Nd and Sm with rela- 
tively small octahedral distortion. In addition, as R goes 
from La, Pr to Nd, Tn gradually decreases with decreas- 
ing spin-exchange constant. In SmTiOs, the TiOe octa- 
hedra are moderately distorted so that the competition 
between the crystal field from the O ions and that from 
the Sm ions comes about. In SmTiOs, the O crystal field 
works as increasing the xy-orhital occupation in all Ti 
sites. As a result, the AFM(G) spin coupling becomes to 
have a two-dimensional anisotropy. The depression of Tn 
of SmTiOs may be attributed to this anisotropy. 

In YTiOs and GdTiOs with large GdFeOs-type 
distortion, the electronic structure is dominantly de- 
termined by the crystal field from the O ions with the 
largely distorted TiOe octahedra. These compounds 
show a FM ground state with an orbital ordering in the 
strong d-type Jahn- Teller distortion. In the previous 
studies, ^'^'-^^^ we have shown that the FM coupling 
in the vicinity of the AFM-FM phase boundary has 
a strong two-dimensional character, resulting in the 
depressed Tq. The puzzling second-order like behavior of 
the AFM-FM transition can be understood from these 
depressions of Tn and Tc- In the previous studies,^''' 
we have also pointed out that although there is no 
stoichiometric system between SmTiOs (AFM(G), 
Ref. 70) and GdTiOs (FM), the AFM(A) state might 
be observed in the solid-solution systems such as 
Lai_xYa;TiOs and Gdi_a;Sma;TiOs sandwiched by the 



Full Paper 



J. Phys. Soc. Jpn. 



Full Paper 



Masahito Mochizuki and Masatoshi Imada 



17 




Ti-O-Ti bond angle 

< 

GdFe03-type distortion 



Fig. 15. Schematic magnetic phase diagram for RTiOs and solid- 
solution systems Lai-j^Yj^TiOs in the plane of temperature and 
the GdFeOs-type distortion. At R =La, Pr and Nd with rela- 
tively small GdFeOs-type distortions, the crystal field (OF) from 
R cations dominantly determines the electronic structures, while 
in significantly distorted GdTiOa and YTiOa, the O crystal field 
dominates. In SmTiOa, both crystal fields strongly compete. 

AFM(G) and FM phases. 

[5]: We finally discuss the experimentally observed re- 
duction of the ordered moment in LaTiOa. The value of 
the ordered moment is about 0.45 /xb, which is reduced 
from spin-only moment of 1 /xb- Since a neutron scat- 
tering study showed that the LS interaction is irrelevant 
to this reduction, the origin has been a puzzling issue. 
Considering that LaTiOs with a small insulating gap of 
-^0.2 eV is located near the metal-insulator phase bound- 
g^j.y^i6,47) have attributed this reduction to the strong 
itinerant fluctuations. Within our effective Hamiltonian, 
the charge fluctuations are completely neglected in the 
limit of the strong Coulomb repulsion. In addition, the 
quantum fluctuations are also neglected in the mean-fleld 
approximation which we have used. Here, let us discuss 
this assignnicint by considering available numerical and 
analytical results on the Hubbard models. 

The Hubbard model is characterized by the ratio U/t. 
The value of U/t for LaTiOs can be estimated as fol- 
lows: The characteristic energy of the nearest neighbor 
Ti t2g transfer t is expressed as Vp^^/A based on the 
second-order expansion with respect to the transfer inte- 
grals between the Ti t2g and O 2p orbitals. The value is 
estimated as 0.33 eV by using the parameter values used 
in our calculation, where A=5.5 eV and T^(i7r=l-33 eV. 
Therefore, the ratio U /t takes approximately 12 for the 
multiplet-averaged Coulomb interaction [/ of 4 eV. 

The equal time magnetic structure factor, which has 
the form; 

^(Q) = ^Ee^^^^'"^'^(KT -ria)K-T -n,-4)), 

(5.47) 

with Q being an AF wave vector, is related to the AF 



ordered moment m through the formula: = -|- 

,f(N). This quantity for the 3D Hubbard model has been 
calculated in several QMC simulations in small-sized cu- 
bic lattices by using the grand canonical method. ^^'^^^ 
However, because of the limitation of the cluster size, 
there exists considerable ambiguity in the extrapolation 
to the thermodynamic limit so that a reliable estimate 
of S{Q) and that of the ordered moment are difficult in 
available calculations. 

At this stage, we discuss these quantities based on 
some approximate results on the 3D Hubbard model 
or numerical results on other models. Then, we present 
some possible mechanisms of the reduction of the ordered 
moment. In particular, we discuss by focusing on effects 
of (1) charge fluctuations, (2) orbital degrees of freedom, 
and (3) imperfect nesting of the Fermi surface. 

We first consider effects of the charge fluctuations. An 
approach based on the Feynman's variational principle 
has been applied to the 3D Hubbard model. This study 
shows that at U/t=12, a relatively large charge fluctu- 
ation remains, which is approximately 40% of that at 
U/t=0. Owing to this charge fluctuation, the ordered 
moment at U/t=\2 is reduced to about 80% of the full 
moment. Here, we note that since this approach neglects 
quantum fluctuations, the ordered moment is saturated 
to 1 /iB in the limit of infinite U (namely, in the limit 
of vanishing charge fluctuation). In the spin- wave the- 
ory on the 3D Heisenberg model, the ordered moment in 
the U ^ oo limit diminishes to ^0.85 /ib owing to the 
quantum fluctuations. Therefore, we naively expect that 
the ordered moment can diminish to ^0.7 ^b when both 
quantum and charge fluctuations are properly taken into 
account. Although the charge fluctuations significantly 
reduce the ordered moment, this value is still a little 
large relative to the observed ordered moment of 0.45 
/Kb, indicating a necessity of additional mechanisms. 

The orbital degrees of freedom can further reduce the 
ordered moment. In the orbitally degenerate systems, the 
AFM spin-exchange interaction is considerably reduced 
as compared with a single band case because of the FM 
contribution from electron transfers between neighbor- 
ing occupied and unoccupied orbitals and the Hund's 
rule coupling. Therefore, the spin correlation in such sys- 
tems is rather weak relative to the nondegenerate sys- 
tems. The resultant relative increase of the charge fluc- 
tuations may play a role on the reduction of the ordered 
moment. Actually, the reduction in the orbitally degen- 
erate system has been studied by Oles.''^) In this work, 
the Hubbard model with flvefold degenerate orbitals and 
isotropic transfers is studied by a Gutzwiller-type varia- 
tional approach, and the ordered moment is calculated as 
a function of {U+6j)/W for various values of j/U where 
j and W are the Coulomb exchange interaction and the 
bandwidth, respectively. If we assume the band width of 
the 3D single-band Hubbard model {W=12t), the val- 
ues of {U+&j)/W and j/U for LaTiOs are estimated 
as ^1.75 and ~0.1, respectively, where U=A eV, j=0.46 
eV and t=0.33 eV. The result shows that for these pa- 
rameters, the ordered moment takes a considerably small 
value of -^0. 45-0. 5 yUe- In LaTiOs, the crystal field from 
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the R cations works as lifting the degeneracy of the t2g 
orbitals at Ti sites. In this sense, this compound may in- 
terpolate between a single-band limit and a completely 
degenerate limit, and the Ti t2g orbitals in LaTiOs are 
also expected to play a role on the reduction of ordered 
moment. 

In addition, we note that when we introduce the next- 
nearest neighbor transfer t' in addition to the nearest 
neighbor transfer t, the AF ordered moment is strongly 
suppressed since t' destroys the nesting of the Fermi sur- 
face. Actually, the reduction due to t' is clearly shown in 
the recent PIRG study on the 2D Hubbard model with 
t and When a small t' of 0.2t is introduced, the 

shape of the Fermi surface is strongly affected, and the 
nested structure is significantly destroyed. As a result, in 
the small U/t region, the ground state is paramagnetic 
and the AFM long-range order appears only in the re- 
gion of U/t >3.45, while for the Hubbard model with 
only t, the AFM long-range order is realized at any posi- 
tive value of U/t due to the perfect nesting. The ordered 
moment m for U/t=4 is estimated as ~0.16, which is con- 
siderably small relative to rn ~0.G for the 2D Hciscnbcrg 
model. More specifically, this value is also much smaller 
than m ~0.35 for the 2D Hubbard model with only t 
at U /t=A7^^ This indicates that the charge fluctuations 
strongly reduce the AF ordered moment when the nest- 
ing of the Fermi surface is destroyed. We expect that 
the GdFeOa-type distortion necessarily reduces the AF 
ordered moment through the destruction of the nested 
Fermi surface. Although the next-nearest neighbor trans- 
fers and the GdFeOa-type distortion are relatively small 
in LaTiOs, deviations from the nested structure together 
with the strong charge fluctuations can reduce the or- 
dered moment. 

On the basis of these discussions, we conclude that 
the reduced ordered moment is qualitatively understood 
from cooperative effects of these mechanisms such as 
charge fluctuations, orbital degrees of freedom and the 
Fermi surface with imperfect nesting. The quantitative 
and microscopic study on the reduction is left for the 
future study. 

6. Summary 

In summary, we have shown that a very simple mech- 
anism of the La crystal field resulted from displacements 
of the La ions caused by the well-known GdFcOs-type 
distortion leads to a consistent picture to dissolve the 
puzzling properties of the AFM(G) state in LaTiOs. We 
have constructed the crystal field Hamiltonians by us- 
ing the experimentally measured position parameters. 
On the basis of these Hamiltonians, we have shown that 
the displaced La ions in the GdFeOs-typc structure gen- 
erates a crystal field with nearly trigonal {Dzd) symme- 
try. Then the threefold degenerate cubic-t2g levels split 
into three isolated levels. We note that the GdFeOa-type 
distortion exists irrespective of the 3d electronic state^^^ 
in contrast with the Jahn- Teller distortion. The lift of 
the Ti Zd degeneracy occurs just as a consequence of the 
crystal field from the La ions. In this sense, the driv- 
ing mechanism proposed here is not a Jahn- Teller but a 
simple and new mechanism, which is inherent and uni- 
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versal with GdFeOs-typc distortion. The energies and 
the spin-exchange constant calculated by the effective 
Hamiltonian show that the lowest-orbital occupations in 
this crystal field stabilize the AFM(G) state, and well 
explains the isotropic spin-wave spectrum with consid- 
erably small spin gap. Our theory have also succeeded 
for the first time in explaining a number of experimental 
findings for LaTiOs in a unified way. The experimental 
results consistently accounted for by our theory include 
neutron scattering, NMR, resonant x-ray scattering and 
x-ray diffraction results. 

In addition, we have also studied the effects of the R 
crystal field in i^TiOs with R being Pr, Nd and Sm with 
relatively small octahedral distortions. We have also re- 
produced the AFM(G) state in these compounds as well 
as the gradual decrease of the spin-exchange constant 
with decreasing size of the R ion, which is consistent with 
the decrease of Tn in experiments. Further, with increas- 
ing GdFeOa-type distortion, the orbital state in the R 
crystal field becomes to have the same symmetry as that 
in the ferromagnetic compounds of YTiOs and GdTiOs 
with large Jahn- Teller distortion, which is in agreement 
with resonant x-ray scattering results. In SmTiOs with 
moderate octahedral distortion, the R crystal field and 
the O crystal field compete strongly, which causes a two- 
dimensional anisotropy in the orbital-spin structure. The 
depression of Tn at R =Sm has been attributed to this 
anisotropy. With decreasing size of the R ion, while the 
R crystal field dominantly determines the orbital-spin 
structure at R =La, Pr and Nd, the O crystal field be- 
comes to dominate over the R crystal fleld at R =Gd 
and Y through the severe competition at R =Sm. As 
observed by several experiments, although the GdFeOs- 
type distortion is rather large, the orbital structures in 
GdTiOs and YTiOs with strong Jahn- Teller distortion 
are well characterized by the d-type Jahn- Teller distor- 
tion, which leads to a FM ground state. 

Recently, it was theoretically predicted that spin liquid 
states with gapless excitation spectra are stabilized in the 
Mott insulators when the geometrical frustration effect is 
large near the Mott transition. As compared with 
the spin degrees of freedom, the orbital degrees of free- 
dom strongly couple to lattice distortions and phonons. 
In general, this difference makes it difficult to stabilize 
the orbital liquid. The present compounds provide typi- 
cal examples of this tendency. 

Through this study, it has been clarifled that in the 
perovskite-type transition-metal oxides, a crystal field 
from the R ions caused by the GdFeOs-type structure 
lifts the 2>d orbital degeneracy. This mechanism of lifting 
the orbital degeneracy competes with the usual Jahn- 
Tellcr mechanism, and eventually plays an important 
role in controlling orbital-spin structures. This mecha- 
nism has general importance since the GdFeOs-type dis- 
tortion is a universal phenomenon, which is seen in a 
large number of perovskite-type compounds. This mech- 
anism may also play important roles on the electronic 
structures in other perovskite compounds. 
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